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The dephasing time of disordered two-dimensional electron gas in modulated magnetic field 
H = (0, 0, H/ cosh 2 ((a; — xo)/S)) is studied. In the weak inhomogeneity limit where 8 is much 
larger than the linear size of the sample, rT" 1 is proportional to H. In the strong inhomogeneity 

limit, it is shown that the dependence is quadratic, tT/ 1 — D (^) H 2 S 2 . In the intermediate 
regime, a crossover between these two limits occurs at H c = ||<5 -2 . It is demonstrated that the 
origin of the dependence of on H lies in the nature of corresponding single particle motion. A 
semiclassical Monte Carlo algorithm is developed to study the dephasing time, which is of qualitative 
nature but efficient in uncovering the dependence of r^> on H for arbitrarily complicated magnetic 
field modulation. Computer simulations support analytical results. The crossover from linear to 
quadratic dependence is then generalized to situation with magnetic field modulated periodically in 
one direction with zero mean, and it is argued that this crossover can be expected for a large class 
of modulated magnetic fields. 



PACS numbers: 73.23.-b, 73.20.Fz, 73.63.-b 
I. INTRODUCTION 

Dephasing is one of the key elements in the study 
of quantum coherent phenomena in mesoscopic systems. 
Coupling with environment suppresses the quantum in- 
terference of electrons. The phase breaking time, or de- 
phasing time is the characteristic time beyond which the 
phase coherence is lost. Although a static magnetic field 
does not destroy all quantum effects, it may introduce 
a cutoff to the interference effect. Dephasing of disor- 
dered electron gas by coupling with a uniform magnetic 
field has been studiedtrcl in the early stage of weak- 
localization theory. In the absence of spin-orbit scatter- 
ing, the magnetic field suppresses the weak-localization 
effect and leads to a positive magnetoconductance, which 
has been observed (for a review, see [4]). The dephasing 
rate due to coupling with uniform field turns out to be 
proportional to the field amplitude. These results can 
also be established ipy qualitative considerations accord- 
ing to KhmclnitskiU. n 

In recent years there has been increasing interests 
in hybrid semiconductor systems both for fundamen- 
tal understanding and for the potentiality of enhanc- 
ing the functionality of the devices. The disposition of 
superconductingJ or magneticu microstructures on the 
surface of heterostructure with a two-dimensional elec- 
tron gas(2DEG) may produce inhomogeneous magnetic 
field which influences electron motion locally. Some in- 
teresting consequences of the modulated magnetic fields 
have already been reported. It was showru that the ef- 
fect of electron-electron Umklapp scattering can be ob- 
served in a 2DEG at a GaAs/Gai-^ALjAs interface by 
imposing a spatially alternating magnetic field normal to 
the 2DEG plane. Kubrak et aliS fabricated a few dif- 
ferent types of hybrid ferromagnetic-semiconductor de- 
vices, which allow them to study how these different 
modulated magnetic fields influence the transport prop- 



erties of 2DEG. And magnetoresistance oscillations due 
to internal Landau band structure of a 2D electron sys- 
tem in periodic, magnetic field are observed^. Theo- 
retically, PengE3 calculated the transport properties in 
a parabolic channel exposed to a periodically modu- 
lated magnetic field. Gumbs and ZhangtU developed a 
magneto-transport theory for the magneto conductivity 
of a square lattice in a periodically modulated magnetic 
field, and predicted some anomalies due to cptnmensu- 
rability effects. Recently Matulis and Peetersia studied 
the semiclassical magnetoresistance in weakly modulated 
magnetic fields. They considered the case when the field 
is periodic in one dimension and with zero mean. In the 
limit of small magnetic field amplitude, it is shown that 
the contribution of the magnetic modulation to magne- 
toresistance increases as J3 3 / 2 in the diffusive limit, while 
increases linearly in B in the ballistic limit. 

The magnetoresistance in the diffusive limit is well de- 
scribed by standard weak-localization theory. It arises 
from the suppression of magnetic field on the cooperon 
propagator which represents the interference of time- 
reversal trajectories. Taking the magnetic fi$ld into con- 
sideration, the cooperon propagator satisfies^, 
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which may be viewed as an imaginary-time Schrodinger 
equation with parameters properly replaced. For the uni- 
form field situation, the solution of the corresponding 
Schrodinger equation is the well-kn^wi Landau levels. 
The dephasing time was found to be 
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where the superscript u represents result for uniform 
field. 

In this paper we first study the dephasing time of dis- 
ordered 2DEG due to coupling with the following modu- 
lated magnetic field 



H = (0, 0, H/ cosh 2 ((a; — xq)/8)). 



(3) 



The corresponding-Schrodinger equation has been solved 
exactly by Hudaklij, where it was shown that a proper 
transformation converts the eouation into a form solved 
earlier by Morsh and Feshbachta. Making use of this solu- 
tion, we construct the cooperon propagator and calculate 
the dephasing time of a disordered 2DEG in this mod- 
ulated field. We find in the weak inhomogeneity limit, 
S >> L, where L is the linear size of the sample, t^ 1 is 
proportional to H, as in the uniform field case. In the 
strong inhomogeneity limit, it is shown that 



1 = £>(-£-) 2 h 2 s 2 . 

T0 \TicJ 



(4) 



In the intermediate regime, a crossover between these two 
limits is expected to occur at H c = §f<5~ 2 - The dephas- 
ing rate dependence on the magnetic field amplitude is 
shown to be related to the nature of corresponding single 
particle motion. Bound states lead to linear dependence, 
while nearly free motion results in quadratic. A semiclas- 
sical Monte Carlo algorithm is developed to study the de- 
phasing time, which is of qualitative nature but efficient 
in uncovering the dependence of on H for arbitrarily 
complicated magnetic field modulation. Computer simu- 
lations support analytical results. The considerations are 
then generalized to situation where the modulated mag- 
netic field is periodic in one direction with zero mean, 
and it is argued that this crossover between linear and 
quadratic dependence can be expected for a large class 
of modulated magnetic fields. 

Before going to quantitative calculations_let's see how 
qualitative considerations of Khmclnitskiitj can predict 
Eq. (Q) . Let the magnetic field be nonzero only in a stripe 
of width S. Consider a large loop of area Dt, we are in- 
terested in the flux piercing it. When 6 is small, the 
effective area where the field is nonzero is S\/Dt. Re- 
quiring the phase change due to this loop to be of order 
1, j^H5yJ Dt^ r~j 1, one immediately obtains Eq.([|). In 
subsequent quantitative calculations, we will see that un- 
derstanding of the origin of this dependence allows us to 
generalize the result to more complicated situations. 

This paper is organized as follows. In section II, the 
solution of the Schrodinger equation by Hudakll3is briefly 
outlined. The construction of cooperon propagator and 
calculation of dephasing time is presented in Section III. 
Section IV contains a description of the numerical algo- 
rithm and the simulation results. The generalization to 
situation with modulated magnetic field periodic in one 
direction is presented in Sec.V. The conclusion is given 
in Section VI. 



II. SOLUTION OF THE SCHRODINGER 
EQUATION IN H= (0, 0, #/ cosh 2 ((V - X )/S)) 

For completeness, the solution of the Schroodinger 
equation in the modulated magnetic field H = 
(0,0, H/ cosh 2 ((a; — Xq) /S)) is briefly outlined in this sec- 
tion. Interested readers are referred to the originaLarticle 
of Hudaklij and the book of Morsh and FeshbachtJ. 

Under the Landau gauge 



A = (0,-ffatanh((x-zo)/<5),0), 
the Schrodinger equation can be written as 



(5) 
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where we have separated variables 



ip(x,y) 



X(x). 



(6) 
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This Schrodinger equation describes the motion of a 
particle in the potential V(x), 



V(x) 



h z 
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where $o = hc/e is the flux quantum, <f> = HS 2 is a 
measure of the magnetic flux by the external field. In- 
troducing 

F = 2tt — , 

P = -pyS/H, 

z = (x- x )/S, 
the potential may be written as 



V(z) 



h 2 F 2 (P 
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For py ^ 0, this potential is an asymmetric well, with, 
different limiting value of y(+oo) and V{—o6). HudakO 
observed that the following transformation 



z/cosh 2 2/i = P 2 + F 2 , 
i/sinh 2 2fj, = 2FP, 



(11) 



can convert the Schrodinger equation into a form that 
was solved by Morsh and FeshbachE£l earlier, 



d 2 X (z) 
dz 2 



+ (e — v cosh 2/i — v sinh 2fi tanh z 
cosh 2 /i \ . 

+ v —- s a) x (z) = o, (12) 

cosh zj 
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where e — 2mS 2 E/h 2 . 

One may distinguish situations between |P| < F and 
|P| > P. In the former case, there is a discrete part in 
the energy spectrum as well as a continuous one, while, 
the latter leads to only continuous part. The solution^ 
is as follows. If 

\P\<F(F-\P\) 2 , (13) 

then the energy spectrum for the motion in the re- 
direction contains a discrete part given by 




with n = 0, 1, [n max ], and 

n max = \jF* + \ -\- y/\P\F. (15) 

The corresponding eigenfunction is 

Xn{z) = iV„exp(-a„z)(e~ z + e z )~ b " x P(-n, 
2^F 2 + 1/4 - n, a n + b n + 1, e~" HfT* + e z )), (16) 

where 



x£ Py (z) = N ktPy exp(- a z)(e~ 2 + e*)- b x 
F(b - 7 + \,b + 7 + \,K+ + 1, e-*/(e-* + e z )), (20) 

with 



a = (k+ + ik)/2, b=(k+-ik)/2, 



l=\jF 2 + -, k + = ^/4F\P\-k 2 . 

This wavefunction vanishes exponentially for x — > oo 
when p y < 0, and bounded from above when x — > — oo. 
For p y > 0, the opposite is true. The inequality ( [l8| ) is 
equivalent to 

< k 2 < A\P\F. (21) 



If the parameters are such that the inequality (^1J) is 
replaced by 

k 2 > A\P\F, (22) 



\P\F 



F 2 + \-{n+\) 



Here F(a, f3, 7, S) denotes the hypergeometric function, 
and N n is a normalization constant. 
Thus there is discrete part if 



(17) 



The continuous spectrum may be divided into two 
parts. For 



(F-|P|) 2 <e<(P+|P|) 2 , 
the energy is given by 



(18) 



(19) 



where k is the momentum in x direction, p(F) = Fh/6. 
The corresponding wavefunction is 



then there is another type of continuous energy spec- 
trum. The wavefunctions are nonvanishing but bounded 
as \x\ — > 00. They describe overbarrier motion. The 
energy is still given by Eq.(19), but the corresponding 
wavefunctions are 



Xk, Py ( z ) 



= N k ^ Py exp(iz(fc 4 
xF((~i(k + k+) - 



-k)/2){e- z + e/' 
l-2 7 )/2,H(fc 



. k + k + 



+l + 2~f)/2,l-ik+,e- z /(e 



K) 
+ e z )), (23) 



where 



Vk 2 



AF\P\ 



The above results are valid for |P| < P. When 
|P| > P, an analytical continuation can be performedlij. 
And it was found that there is no discrete part in the 
spectrum, the wavefunctions and energy are of the same 
form as that for |P| < P, but now valid for this region of 
parameters. 
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III. DEPHASING TIME 

In the weak inhomogeneity limit i5 > L where L is 
the linear size of the sample, F> 1 for not very weak 
field. Then the potential well (given by Eq.(||)) is deep 
enough to host many discrete levels, which are reminis- 
cence of the Landau levels. Recall that low energy states 
dominate in the cooperon propagator, the continuous 
part (with energies higher than the barrier height) gives 
negligible contribution, thus only the discrete part of the 
spectrum is important in the cooperon. It can be shown 
that in this limit, inhomogeneity brings a correction of 
order 0{y) to the usual uniform field weak- localization 
magnetoresistance. And the dephasing time has the same 
form with the situation of a uniform field. The opposite 
limit, F <C 1, which we will focus on, is more interesting. 

Assuming the field described by Eq.(||) can be real- 
ized in experiment, let's show that for realistic parame- 
ters it is possible to have the discrete part of the spec- 
trum absent in the strong inhomogeneity limit. For 
S ~ lOOnm, F ~ ^f- ~ IQH if H is in Tesla, thus 
up to H ~ lOOGauss, one can take F as a small quan- 
tity. Then we check if the discrete part of the spec- 
trum exists in this limit. The criterion for its existence 
Pvl ^ p3 ^ m-3 



is p = E|£ < f 3 < 1(T 3 . Since p y > one has 
> j-. For a system with La, < 10 4 nm, this 



P 



10 . Thus the inequality for 



„ ives P _ Pv£ > io>" 

gives r h 10 4 nm 

the existence of discrete levels does not hold in this situ- 
ation. One therefore concludes there is no discrete part 
in the spectrum. In this limit, 



E : 



[pl + {\Pv\-Ff 



2mS 2 

and the cooperon propagator is 



(24) 



C(r,r';w) 
6 



dp x dp y 



iwT+^\pl + {\ Pv \-Ff 



(25) 



where the integral should be done under the constraint 



Dt 



— [pl + (\p y \-Ff] «1, 



P., 



(\ Py \-Ff« q 2 , 



(26) 



7^ (27) 



Dt 



as it is a condition for the perturbation theory. 

Before proceeding, we need to discuss the boun dary 
condition. If we use the wavefunctions Eq.(^TJ) and d23| ) 
in their present form, then the cooperon propagator con- 
structed naively in Eq.(p5f) docs not satisfy the zero am- 
plitude boundary condition 



C(±L,y) = C(x,±L) = 0. 



(28) 



The boundary condition in y direction can be taken care 
of in the same way as in the uniform field case, where 
the plane wave wavefunctions are replaced by their lin- 
ear combinations, and p y takes only some discrete allowed 
values. The boundary condition in the x direction can 
be dealt with in the same manner. We note that from 
Eq. (|l9|) , the eigenstates with k and — k are degenerate, 
thus a linear combination of them is also an eigenstate. 
The new eigenstate can be constructed as 



<:*M = <: B Py +Px%l Py , (29) 



which are required to satisfy the boundary condition 

(£)=<£,(-£) =o. (30) 

As in the case of uniform field, only some discrete p x will 
be allowed. Thus the continuous part of the spectrum 
becomes again discrete in this boundary condition. The 
integration over momentum will be replaced by summa- 
tion over these discrete allowed values. We note that 
the dispersion relation is not altered in the construction 
Eq.©. 

Introducing the dimensionless momenta q x = p x /F, 
q y = Py/F, q' = q/F, P' = P/F, and going to time 
domain, we find 



C(r,r';t,t') 
h 2 F 2 



5 2 i 



E 



^L, gtl (x,y)^q x , qv (x',y') x 



x exp 



ql+(\qy\-i) 2 <q' 
DF 2 



S 2 



<zI + (M- 1) 2 ] (*'-*)} (3i) 




E 



i>f* qy (x,y)ipf xtqy (x',y') 



ql>MP'\ 



x exp 



^ + (|9«|-1) 2 <9' 2 



DF 2 



[<& + {\%\-lf] (t'-t) 



We see that the magnetic field results in a characteristic 
time scale beyond which the cooperoon propagator or the 
interference effect is no longer important. Therefore 



1 

T 4> 



DF 2 



(32) 



which leads to Eq.(@). 

Thus the dephasing rate depends quadratically on the 
field amplitude in the strong inhomogeneity limit. In this 
calculation, we see that the dependence of on H is a 
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result of the nature of corresponding single particle mo- 
tion. From Eq.(|8|), when F is large, the potential well 
is deep enough to accommodate many bound states with 
discrete levels(which are reminiscence of Landau levels). 
Dominant contribution to the cooperon is from these low 
lying levels which results in linear dependence. In con- 
trast, if F is small, the potential well can be so shallow 
that discrete level does not appear, and the particle ex- 
ecutes overbarrier motion which is nearly free with con- 
tinuous spectrum. When this part of spectrum is domi- 
nant in the cooperon, the dephasing rate is a quadratic 
function of the field amplitude. The quadratic depen- 
dence implies that the electrons are more slowly dephased 
compared with in the uniform field case, where the de- 
pendence is linear. Physically this is because for a given 
diffusion time, electrons in this regime can visit a larger 
area than the constrained motion. Since the magnetic 
field is nonzero only in a limited region, the phase change 
accumulated during this time is smaller than that of the 
uniform case. So the electron has to wander for a longer 
time before it gets dephased. Therefore the constrained 
motion of single particle leads to linear dephasing rate 
dependence on field amplitude, and nearly free overbar- 
rier motion to quadratic. 

Eq.(||) is obtained in the limit of small F, which may 
be realized with tiny S and small H. More realistic 
for experimental observation is situation with moderate 
6. In this case, one expects a crossover from the lin- 
ear(uniform limit) to quadratic (inhomogeneity limit) de- 
pendence. This point can be illustrated in the following 
way. Consider the electron which has diffused for a time 
t. If Dt < 5 2 , the inhomogeneity is not noticed by the 
electron. If during this time the electron has already 
been dephased, r > r?, where is the dephasing time 
in the uniform field(given by Eq.(j^)), then the inhomo- 
geneity is not important at all. When r is beyond this 
time scale, the inhomogeneity effect enters. Thus, for a 
given S, there is a crossover field implied by 

(33) 



which yields 



He. = 



he 



4eP- (34) 
This is equivalent to saying that the inhomogeneity effect 
enters when F *~ 1. 

On the other hand, if we can adjust the width S, then 
for a given H, there is a crossover length scale S c , 



5,= 



he 

4efF' 



(35) 



which is the magnetic length for a uniform field. 

In passing on, let's mention that from Eq.(|3l]), one 
can see that the magnetoconductance is proportional to 
H 2 8 2 in this small F limit. In the large F limit, a log- 
arithmic dependence is expected as in the uniform field 
case. These can also ripe established by qualitative argu- 
ment of Khmclnitskiila. 



IV. MONTE CARLO SIMULATION 



A Monte Carlo algorithm has been developed to sim- 
ulate the dephasing process. In order to be consistent 
with that implied by Eq. (1), the simulations will be 
semiclassical in nature. Trajectories will be used and the 
only quantum mechanical effect will be in the phases. 
In this approach, a particle performs random walk in a 
square lattice. The value of the perpendicular field is as- 
signed to a dual lattice. We trace all the closed loops that 
are formed. Once a loop is formed, all the inner points 
are picked up(for technical details, see [17]), and we cal- 
culate and record the phase change fa due to this loop. 
Then the trajectory of this loop is erased, and the parti- 
cle continues the random walk. The phase accumulated 
this way is 



(36) 



whose average is zero since the loop has equal possibility 
to be clockwise and counter-clockwise. Then 



(37) 



The random walk stops as soon as 



b 2 ) > 1, 



(38) 



and we specify the time(total number of random walk 
steps) that has been spent to reach this as the dephasing 
time. An ensemble of such walks are performed and the 
dephasing time is averaged. In the simulations, the mean 
free path I and the flux quantum <3>o have been set to 1, 
so that 



(39) 



where j runs over all the inner points in the dual lattice 
of loop i. 

Several features make this simulation of qualitative na- 
ture. First, the criterion Eq.([38|) is a qualitative one. 
Secondly, in the simulation, the impurities are assumed 
to be on regular lattice sites, with the mean free path 
as lattice constant. Thirdly, sometimes we choose the 
field not necessarily the original one, but qualitatively 
the same. However, essential physics is not lost despite 
of its qualitative nature, and this simulation is very ef- 
ficient in uncovering dependence of the dephasing time 
on the field amplitude, which is particularly appropriate 
for situations where the field modulation is complicated 
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so that it is difficult to make progress with analytical 
approaches. 

As a check of the algorithm, we have simulated the uni- 
form field case. Excellent linearity in the plot of dephas- 
ing rate against H is observed, in agreement with Eq. (^|) . 
For comparison with the modulated field we studied in 
the previous section, the simulation is performed with 
the following field 



H(x) 



H, \x\ < 5, 
0, otherwise, 



(40) 



which is qualitatively the same as Eq.(l). The results 
plotted in Fig. ([!]) show a clear linear dependence for large 
<5, and a crossover from linear to quadratic dependence 
for moderate S. The boundary is set at L — 5000/ in this 
simulation. Trajectories which touch the boundary are 
excluded. The results are obtained by averaging over 10 4 
random walks. When the inhomogeneity effect enters, 
fluctuations become significant in the result. This is be- 
cause there are extreme trajectories that wander for a 
long time in the zero field region. However, by recording 
not only the average but also the standard deviation, we 
are able to get an estimation 2.0 ±0.1 for the exponent 
in the inhomogeneous limit, in agreement with Eq.(^|). 
Thus the numerics support analytical results. 



□ 5 = 100Z 
• 8 = 2001 
o.ooi - o 5 = 1000Z 




H 

FIG. 1. Simulation results for the dephasing time in the 
modulated mag netic field Eq.@. The data are obtained 
by averaging over 10 4 random walks. The three curves are 
for S = 1000/, 200/, 100/, respectively. The boundary is 
L = 5000/. The straight line is a fit to linear dependence, 
dotted line that of quadratic dependence. 



V. GENERALIZATION TO MODULATED 
MAGNETIC FIELD PERIODIC IN ONE 
DIRECTION WITH ZERO MEAN 

The qualitative feature of the results, linear depen- 
dence in weak inhomogeneity limit and quadratic in 



strong inhomogeneity limit, may be more general. Con- 
sider the situation with magnetic field modulated peri- 
odically in one direction with zero mean. While the weak 
inhomogeneity limit is easy to understand, we shall fo- 
cus on the opposite limit to discuss possible quadratic 
dependence of 1/r^ on H. 

Assume the magnetic-field profile is such that near- 
est neighbor stripes of width 5 have the same magnitude 
but opposite sign. Consider the phase change of a closed 
loop due to this magnetic field. Since the field is periodic 
with zero mean, what is important for the magnetic flux 
is the ratio of linear size of the loop with the width of one 
stripe, yjDr^/S. If this ratio is an even integer, then the 
net flux is zero. If it is an odd integer, then the flux is 
equal to that piercing through a single stripe. In general 
this ratio may fluctuates, however it is clear that the av- 
erage flux is proportional to y/Dr^Hd. Demanding the 
phase change due to this flux to be 1, one obtains Eq.([|), 
with probably an additional numerical coefficient. 

For a general magnetic field periodic in one di- 
rection, the corresponding Schrodinger equation is 
often difficult to solve analytically. For example, 
for H = (0, 0, H sin(x/S)), one can choose A = 
(0, — H5 cos(x/S), 0). Then the resulting Schrodinger 
equation is a Whittaker-Hill equation whose solution may 
be reduced to three-term recurrence relations!!^, but an 
explicit analytical result for the spectrum is not known. 
However some insight can be gained by analyzing the 
structure of energy spectrum. In general, if the magnetic 
field is (0, 0, H(x)), in the Landau gauge A = (0, A(x), 0). 
Separating variables in the wavefunction, the electron 
motion in x direction is described by a Schrodinger equa- 
tion with potential 



V(x,p v ) = p v - -A(x) 



(41) 



which is also a periodic function. For given p y , this poten- 
tial is bounded from above as that described by Eq.(||). 
The profile of this potential is a series of potential wells 
joining with each other. The energy spectra of electrons 
in a family of such modulated-piagnetic fields are calcu- 
lated by Ibrahim and Peeteralj. For electron energy less 
than the barrier height, the wavefunctions in neighboring 
wells overlap and spread. So the discrete levels in single 
wells now form minibands due to the periodic structure. 
When H and S are large, the barrier can be quite high, 
and states in these minibands dominate in the cooperon 
propagator. According to the experience in dealing with 
the previous case, the dephasing rate is expected to be a 
linear function of H . When H and 8 are small, however, 
the potential well is not deep enough to support these 
minibands, and dominant contribution is expected from 
nearly free overbarrier motion, which leads to a quadratic 
dependence. 

To examine these considerations, computer simulations 
are performed for the following field 



G 



H(x) 



(42) 



H, 2nS < x < (2n + 1)5, 
-H, {In - 1)5 < x < 2nS. 

The results are shown in Fig.|| The curves show much 
resemblance to that in Figjlj, and a crossover from lin- 
ear to quadratic dependence is evident. There are also 
some differences which will be noted here. First, for a 
given 5, the crossover field H c is smaller in the field ( [42] ) 
than in ( [40|), which suggests when H is small, electron in 
the field (|42|) is more quickly dephased than in the field 
(40). This may look strange at first sight, since the field 
(42) has zero mean, one may expect the electrons in this 
field are more slowly dephased. However, this result can 
be understood by noticing what enters is the variance of 
the phase instead of the average, and the field ( f42| ) is 
nonzero everywhere in the plane. Another difference is 
in the field ( f42| ) there is broader crossover region. This 
is associated with the fact that, in the periodic situa- 
tion, states with energy lower than the barrier height are 
not really bound states. They share some feature of the 
plane waves according to Bloch's theorem. As discussed 
previously, this feature has the tendency of leading to a 
quadratic dependence. When 5 is moderate, the barrier 
is not very high for many values of p y . Then these states 
can be dominant in the cooperon, resulting in a broad 
crossover region. 



o S = 201 
• S = m 

d8 = 4001 
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FIG. 2. Simulation results for the dephasing time in the 
modulated magnetic field periodic in one direction given 
by Eq.(^). The data are obtained by averaging over an 
ensemble of 10 3 random walks. The three curves are for 
8 — 400Z, 40/, 20Z, respectively. The straight line is a fit to 
linear dependence, the dotted line that of quadratic depen- 
dence. 



that the crossover from a linear to quadratic dependence 
in the dephasing rate on field amplitude can be quite 
general. For situations where the corresponding single 
particle potential is bounded from above, the particle 
motion can exhibit both bound and free nature, leading 
to a crossover from linear to quadratic dependence. 



VI. CONCLUSION 

We have studied the dephasing time of disordered two- 
dimensional electron gas in modulated magnetic field 
H = (0,0, H/ cosh 2 ((x — xq)/S)). It is shown that in 
the weak inhomogeneity limit, 5 >> L, where L is the 
linear size of the sample, r^ 1 is proportional to H . This 
happens when the bound states with discrete spectrum of 
the corresponding Schrodinger equation dominate in the 
cooperon propagator. While in the strong inhomogeneity 
limit, the dependence is quadratic t^ 1 = D (^) H 2 S 2 . 
In this case, the nearly free overbarrier motion gives 
dominant contribution to cooperon. In the intermediate 
regime, a crossover between these two limiting situations 
occurs at H c — ^5~ 2 . A semiclassical Monte Carlo al- 
gorithm has been developed to study the dephasing time, 
which is of qualitative nature but efficient in uncovering 
the dependence of on H for arbitrarily complicated 
magnetic field modulations. Computer simulations sup- 
port analytical results. These considerations are general- 
ized to the situation with magnetic field modulated peri- 
odically in one direction with zero mean, where a similar 
crossover is observed. We believe this crossover between 
linear and quadratic dependence can be expected for a 
large class of modulated magnetic fields where the corre- 
sponding single particle potential is bounded from above, 
so that the motion exhibits both bound and free nature 
depending on parameters. 
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Numerical experiments are also performed for situa- 
tions with some other magnetic field modulations, for 
example, that with circular symmetry, or periodic in 
both directions. Similar crossover behaviors have been 
observed. These simulations as well as the qualitative 
understanding of the origin of the dependences, suggest 
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